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Abstract: The number of galaxies at a given flux as a function of the redshift, z, is derived 
when the z-distance relation is non-standard. In order to compare different models, the same 
formalism is also applied to the standard cosmology. The observed luminosity function for 
galaxies of the zCOSMOS catalog at different redshifts is modelled by a new luminosity 
function for galaxies, which is derived by the truncated beta probability density function. 
Three astronomical tests, which are the photometric maximum as a function of the redshift 
for a fixed flux, the mean value of the redshift for a fixed flux, and the luminosity function 
for galaxies as a function of the redshift, compare the theoretical values of the standard 
and non-standard model with the observed value. The tests are performed on the FORS 
Deep Field (FDF) catalog up to redshift z = l.b and on the zCOSMOS catalog extending 
beyond z = A. These three tests show minimal differences between the standard and the 
non-standard models. 
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1. Introduction 

The linear correlation between the expansion velocity, v, and di, the distance in Mpc, is 


V = Hodi = cz 


( 1 ) 
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where Hq is the Hubble eonstant Hq = lOO/ikms ^ Mpc with h = 1 when h is not speeified, c is the 
veloeity of light, c = 299792.458^ see Mohr et al. [1], and z is the redshift defined as 

^ohs ^em 

^ i’ (2) 

'^em 

with Xobs and Agm denoting respeetively the wavelengths of the observed and emitted lines as determined 
from the lab souree, the so ealled Doppler effeet. This linear relation ean be derived from first 
prineiples, namely general relativity (GR), and is only a low-^ limit of a more general relation in standard 
eosmology. The presenee of a veloeity in the previous equation has pointed the standard eosmology 
towards an expanding universe. The previous equation also eontains a linear relation between distanee 
and redshift, and this has pointed the experts in plasma physies towards an explanation of the redshift 
in terms of the interaetion of light with the eleetrons of the intergalaetie medium (IGM). In terms of 
plasma physios, the expansion veloeity of the universe beoomes a pseudo-veloeity beeause the Universe 
is supposed to be Euelidean and statie. The eonjeoture here adopted is now outlined: the redshift should 
be oonsidered as a speetroseopie measure independent of the adopted eosmology. An example of this 
eonjeeture ean be found at the home page of the Supernova Cosmology Projeet (SCP) ‘All the analyses 
were developed with eosmology hidden.’ The main physieal explanations for the redshift are: a Doppler 
shift, whieh means an expanding universe, a general relativistie effeet, see [2], a plasma effeet, see [3], 
and a tired light effeet as suggested by [4,5]. More details on the various theories whieh explain the 
eosmologieal redshift ean be found in [6]. A point of discussion: The presenee of physieal effeets whieh 
explain the redshift allows speaking of a Euelidean universe in whieh the distanees are evaluated with 
the Pythagorean theorem. In a Euelidean universe, the main parameters are Hq, z, and the distanee, d, 
but the veloeity of expansion is virtual rather than real. The aim of having a Euelidean universe is the 
explanation of the astronomieal variables sueh as the redshift and the absolute magnitude and eount of 
galaxies without GR. 

Coneeming the value of Hq and Dm, we will adopt reeent values as obtained by a mixed model whieh 
uses the eosmie mierowave baekground (CMB) measurements at high redshift and the baryon aeoustie 
oseillation (BAG), see [7], 

/7o = (69.6 ± 0.7) kms“^ Mpc“^ , (3) 

Dm = 0.286 ± 0.008 . (4) 

Hubble’s eonstant is explained in the dynamieal relativistie models beginning with [8,9]. Reeently, 
researeh in the framework of modem theories on an aeeelerating universe has been foeussed on 
measuring eosmologieal parameters sueh as Dm and Ua, see [10,11]. In the last years, the enormous 
progress in astronomieal observations has inereased the available data for galaxies up to 2 ; =3.36, see the 
EORS Deep Eield (EDE) eatalog, whieh is made up of 300 galaxies with known speetroseopie redshift, 
see [12,13]. Another high redshift eatalog is zCOSMOS, whieh is made up of 9697 galaxies up to 2 : =4, 
see [14]. These data demand a new formalism for the number of galaxies as a funetion of the redshift. At 
the same time, Hubble’s law ean be inserted into a more preeise physieal framework. In order to eover 
these questions, Seetion 2 reviews some old and new derivations of Hubble’s law as well the magnitude 
system, and Seetion 3 derives a new relation for the number of galaxies as a funetion of the redshift. 
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Section 5 introduces a relativistic model for the number of galaxies as a function of the redshift and 
Section 6 deals with the luminosity function for galaxies at different redshifts. 


2. Basic Formulae 


The change of frequency of light in a gravitational framework [15], a photo-absorption process 
between the photon and the electron in the intergalactic medium, see Equation (3) in [4], or a plasma 
effect, see Equation (50) in [3], all give 

V = c [exp(ifo ~ 1] • (5) 


We can isolate the distance di in Eq. (1), obtaining the linear relation 


di = ^ . ( 6 ) 

-no 

In standard cosmology, one cannot use the previous approximation when inserting d into Eq. (10) for 
redshifts larger than 0.1. The previous equation needs to be generalised in standard cosmology when 
z > D.l with the use of both the ‘Hubble function’, see Eq. (58), and the luminosity distance, see 
Eq. (62). 

The expression for the distance, d, in the nonlinear Eq. (5) gives the relation 


d 


In {z + l)c 
^0 


A Taylor expansion around z = 0 of this equation gives 


d 


cz 1 1 cz^ 

iro~2iro^3H^ 


(7) 


( 8 ) 


In the limit z —?• 0, the linear distance, di, and the nonlinear, d, are equal. The laboratory results of 
the line shift in dense and hot plasmas can be found in [16-20]. As a first example, the experimental 
verification of the redshift of the spectral line of mercury as due to the surrounding electrons can be 
found in Eigure 1, see also [21]. A second example is given by the Balmer series Ha line emitted by 
laser-produced hydrogen. A linear fit of the data in Table 1 in [23] gives the relation for the redshift of 
the Ha line 

z = 0.0042 + 0.00039A^e when 0.0042 < z < 0.0568 , (9) 


where A’e is the electron number density. This linear relation can be visualized in Eigure 2. The previous 
relation shows a laboratory redshift as a function of a growing density. In astrophysical plasma, we 
are interested in the redshift as a function of the distance when the temperature of the medium is 3K, as 
known from the cosmic microwave background (CMB), and its mean density is extremely low. Therefore 
the previous laboratory experiment is illustrative rather than quantitative. 
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Figure 1. Hgl 435.83 nm line shifts versus the eleetron density, data as extraeted by the 
author from Figure 7 in [22] (empty stars) and linear regime (full line). 



Figure 2. Redshift of the Ha line versus the eleetron density, data as extraeted by the author 
from Table 1 in [23] (empty stars) and linear fit (full line). 
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Figure 3. The blue absolute magnitude Mb computed with the nonlinear Eq. (11) for 
263 galaxies belonging to the PDF catalog versus the well measured redshift. The lower 
theoretical curve as represented by the nonlinear Eq. (11) is the red thick line when 
mL=30.33, which is the maximum apparent magnitude of the catalog, M.q = 5.48 and 
Hq = 69.6 km Mpc“^ (green points). The redshift covers the range [0,4.5] 

2.1. Magnitude system 


The absolute magnitude of a galaxy, M, is connected with the apparent magnitude m through the 
relation 

M = m-51ogio(d)-25 . (10) 

The nonlinear absolute magnitude as a function of the redshift as given by the nonlinear Eq. (7) is 

dn {z + 1) c'' 


M = m — 5 log 


10 




25 


( 11 ) 


The previous formula predicts, from a theoretical point of view, an upper limit on the maximum absolute 
magnitude which can be observed in a catalog of galaxies characterized by a given limiting magnitude. 
The previous curve can be connected with the Malmquist bias, see [24,25], which was originally applied 
to the stars and later on to the galaxies by [26]. We now define the Malmquist bias as the systematic 
distortion in luminosity or absolute magnitude for the effective range of galaxies due to a failure in 
detecting those galaxies with fainter luminosity or high absolute magnitude at large distances. 

Figure 3 shows such a curve, the Malmquist bias, for the PDF catalog. The well measured 
spectroscopic redshift, the blue and red apparent magnitude of the PDF data, can be found at the 
Strasbourg Astronomical Data Centre (CDS). Figure 4 shows the curve of the Malmquist bias for the 
zCOSMOS catalog. A careful examination of Figures 3 and 4 allows to conclude that all the galaxies are 
in the region over the border line of the Malmquist bias and the number of observed galaxies decreases 
with increasing z as theoretically predicted. 

In a Euclidean, non-relativistic, and homogeneous universe, the flux of radiation without attenuation, 
/, expressed in units of where Lq represents the luminosity of the sun, is 


/ = 


L 




( 12 ) 
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Figure 4. The B absolute magnitude M computed with the nonlinear Eq. (11) for 9697 
galaxies belonging to the zCOSMOS catalog versus the redshift. The lower theoretical curve 
as represented by the nonlinear Eq. (11) is the red thick line when 771^=23.2, A4q = 4.08 
and Hq = 69.6 km Mpc“^ (green points). The redshift covers the range [0,1] 

where d represents the nonlinear distance of the galaxy expressed in Mpc, see Eq. (7) 

connecting the absolute magnitude, M, of a galaxy to its luminosity is 

— = 

r ’ 

Lq 

where Mq is the reference magnitude of the sun in the considered bandpass. 

The flux expressed in units of as a function of the apparent magnitude is 

/ = 7.957 X 10® m 

Mpc^ 

and the inverse relation is 

m = Mq - 1.0857 In (0.1256 X 10“®/) . (15) 

Once a band is fixed, we have a reference magnitude of the sun in that band. We give a few examples: 

band u* in SDSS Mq=6.3S, band B in EDE Mq= 5.48 and band B in zCOSMOS Mq=4.08. 


The relation 

(13) 

(14) 


2.2. Tired light 

Assume that the photon loses energy, E, in a way proportional to its energy: 

dE 

—— = —CostE . 
dt 

The coefficient Cost is assumed to be proportional to the averaged number density of the IGM, n. 


(16) 


Cost = an , 


(17) 


and therefore 


dE 

—— = —anE 
dt 


( 18 ) 
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We now replace the energy E = hu where h is Planck’s constant, 

dv 


dt 


= —anv 


and we convert the frequency to the wavelength. 


dX 

— = —andt 
A 


(19) 


( 20 ) 


The speed of light, c, is assumed to be constant, % = c 


dt 


, A a 


nds 


( 21 ) 


where Aq is the original wavelength. On introducing the redshift 

A — Ao 


z = 


An 


we obtain 


ln(l + z) = - < n > d , 
c 


( 22 ) 


(23) 


where < n > is the averaged number density. The Hubble constant can be introduced with the following 
meaning, 

Hq = a < n > , (24) 


obtaining 


ln(l + z) = —-d . 
c 


(25) 


The distance d has the same dependence on the photo-absorption/plasma process as given by Eq. (7) 

In {z + l)c 


d = 


Ha 


(26) 


The observed flux without the absorption of our galaxy is 




(27) 


and the distance modulus for tired light is 


m — M = 25 + 5 


In (10) 


In 


In (1 + z) cA 5 In (1 + z) 
Ho ) ^2 In (10) 


(28) 


A generalization of the concept of tired light can be obtained by inserting a supplementary absorption of 
the light, i.e., Compton scattering, see formula (51) in [3]: 


, i , Ho,, „Ho„ L 1 


(29) 


where /) is a variable parameter which is 1 when only tired light is considered and 3 when the Compton 
scattering is added. Here we have invoked the Compton scattering as a possible source of absorption but 
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the parameter fi ean be eonsidered a regulating parameter of an unknown seattering meehanism. The 
distanee modulus of generalized tired light without galaetie extinction is 


m — M = 25 + 5 


In (10) 


In 


In {1 + z)c 

Wn 


5 /3 In (1 + ^) 

^ 2 In(10) 


(30) 


Hubble’s constant can be extracted from this equation as a function of the distance modulus (m — M): 

f/o = 100000 In (z +kms"^Mpc"^ , (31) 


and, as a practical example, when m — M=43.834 and 2 ; = 0.974, which is the case with SN C-001 in 
the Union 2.1 compilation, we have Hq = 69.6 kms“^ Mpc“^ when = 2.032, which is the same value 
quoted in our Introduction. 

A first comparison can be done with the distance modulus in a plasma environment as given by Eq. (7) 
in [27] without galactic extinction: 


m — M 


In (In (2: + !)) ^ ln(^ + l) 1 m A V 95 

In (10) ^2 In (10) ^ In(lO) V^o/ 


(32) 


see equation (7) in [27]. A second comparison can be done with the historical equation (23) in [28], 
which is the [m, z] relation for the steady-model 


m — M 


In(lO)^ lii(lO) ^ lii(lO) \Ho) 


+ 25 


(33) 


We briefly recall that Sandage [28] used two explosive models, a Friedman model [8,9] and a steady 
state model [29] . 

Figure 5 presents the behavior of the three distance moduli here considered as functions of the redshift. 
We now outline how formula (26) dates back to the year 1929. The original formula for the change of 
frequency of light in a gravitational framework is due to [15]: 

Ai/ lAnGpDL 

— = -+- • (34) 

Here, v is the frequency, G is the Newtonian gravitational constant, p is the density in g/cm^, D is the 
distance after which the perturbing effect begins to fade out, L is the distance, and c is the speed of light. 

We make the change of variables ^ ^ ^ = —dE/h, and u = E/h, 

which means 

dE a , 

—— = — < n > ds . (35) 

E c 

On repeating the previous procedure, we obtain formula (26). 

There is actually a debate on the existence of tired light in laboratory experiments. Here we use tired 
light as a simple theory to model a more complex interaction between light and the intergalactic plasma. 
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Figure 5. Distance modulus for the generalized tired light as given by Eq. (30) (full line) 
with (5=2.1, for a plasma as given by Eq. (32) (dashed line), and for the steady-model as 
given by Eq. (33) (dot-dash-dot-dash line). 

2.3. The luminosity function 

The Seheehter funetion, introdueed by Seheehter [30], provides a useful fit for the luminosity funetion 
(EE) for galaxies: 

«>(E)dL = (£:)(E)Oexp(-T)rfi, (36) 

Here, a sets the slope for low values of L, L* is the eharaeteristie luminosity, and is a normalization. 
The equivalent distribution in absolute magnitude is 

$(M)dM = exp(-10°-^(^*-"))dM, (37) 

where M* is the eharaeteristie magnitude as derived from the data. The sealing with h is M* — 5 logj^g ^ 
and [Mpc~^]. 


3. N-z relation 


This seetion evaluates the number of galaxies as a funetion of the redshift, firstly assuming a linear 
relation and seeondly a nonlinear relation between redshift and distanee. The evaluations are done on a 
sphere of radius r whieh is identified with the ehosen distanee. The main statistieal test is the x^' 



(38) 


where n is the number of bins, Tj is the theoretieal value, Oi is the experimental value represented in 
terms of frequeneies, and cXi is the error eomputed as the square root of Oi. 
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3.1. The linear case 


The joint distribution in z and/for galaxies, see formula (5.133) in [31] or formula (1.104) in [32] or 
formula (1.117) in [33], is 


dNi / C 4. . ^ , 

= 47r(—) , 


dVtdzdf 


(39) 


"'crit 


where dfl, dz and df represent the differential of the solid angle, the redshift, and the flux, respectively, 
and <f) is the Schechter LF. The critical value of 2 ;, Zcru, is 




H^L* 

Anfc'^ 


(40) 


The number of galaxies in 2 : and / as given by formula (39) has a maximum at 2 ; = zi^pos-max, where 

^l,pos—max ^crit \/cr T 2 , (41) 

which can be re-expressed as 


^Lpos—maxi.f') /— n i 

2 0rV/c 

or introducing the two observable variables, Mg^n and M*, 

1.772 X 10-5 V2Ta\/ 10°'^^©"°'^^*i7o 


(42) 


^l,pos—max(.'^') 


0.921 -0.921m 

^TTV e W C 


(43) 


These two formulas, which model the photometric maximum, are not reported in chapter 5 of [31]. 
The theoretical mean redshift of galaxies with flux /, see formula (1.105) in [32], is 

r(3 + a) 


^l,crit 


r(5/2 + a) 


(44) 


3.2. The nonlinear case 


We assume that / = and 


= ^{Horlc) _ Y ^ 


(45) 


where r is the distance; in our case, d is as represented by the nonlinear Eq. (7). The relation between dr 
and dz is 

dr = 


cdz 


{z + l)H^ • 

The joint distribution in z and/for the number of galaxies is 

= h I 

where 5 is the Dirac delta function. 

The evaluations of the integral over luminosity and distances gives 




(46) 


(47) 


dN 


dD.dzdf 


= F{z-f,n) = 


A /n { (ln(z + l))2 

4 (In [z + 1))^ f I 6 ^cri3 TT 


Ho^L* (z + 1) 


(48) 
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The number of galaxies in z and / as given by formula 48) has a maximum dX z = Zpos-max, "where 

^ Lstar Hq — Hq 


Z (f)=e-^^^^ 

'^pos—max\J ) ^ 

or introducing the two observable variables, Msun and M 


■K C^f 


1 , 


(49) 


e-i/16 


^pos—maxi^^') 


^]^q0.4 Msun-0.4 //p _ ^@4 ,r a c2/+128 ir c2/ + io0.4 Msun-0.4 M* Msun-0.i M* ^ 

n c^f — 


(50) 


The total number of galaxies comprised between a minimum value of flux, fmin, and a maximum value 
of flux fmax, for the Schechter LF can be computed through the integral 


dN 

dildz 



Fiz;f,n)df 


(51) 


This integral does not have an analytical expression and we must perform a numerical integration. 
The theoretical mean redshift of galaxies with flux / can be deduced from Eq. (48): 


f,-F{z-f,n)dz 


(52) 


The above integral does not have an analytical expression, and should be numerically evaluated. The 
differences between the formulas of this subsection and the formulas in chapter 5 of [31] are that our 
formalism is built to cover high values of the redshift, against the low values of the redshift of the 
standard approach. Both models are built in the framework of a Euclidean universe. 


4. Astrophysical applications 

We processed two catalogs in order to test the theoretical formulae: the EDE and the zCOSMOS. 
These two catalogs differ for the number of galaxies and parameters of the Schechter luminosity function 
(EE) for galaxies. We now review the three parameters of the Schechter EE: a is fixed, is not relevant 
because we equalize the theoretical maximum in frequencies, and M* is allowed to vary in order to 
match the observed frequencies. The observed mean redshift of galaxies with a flux / or m, (zobs), is 
evaluated by the following algorithm. 

1. A window in apparent magnitude or flux is chosen around m or /. 

2. All the galaxies which fall in the window are selected. 

3. The mean value in redshift of N selected galaxies is (zobs) and the uncertainty in the mean, cr^, is 
ap = s/\fN where s is the standard deviation, see formula (4.14) in [34]. 
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data 


Figure 6. The galaxies of the PDF eatalog are organized in frequeneies versus speetroseopie 
redshift. The redshift eovers the range [0,4] and the histogram’s interval is 0.1. 

4.J. The PDF catalog 

The peneil beam eatalog FDF has a solid angle of ~ 5.6 sq aremin, or T x T around the south 
galaetie pole, and eovers the interval [0,4] in redshift. In partieular, we seleeted the 263 galaxies 
with speetroseopieal redshift and we proeessed the B band whieh has the range in apparent magnitude 
[19.5, 30.2] mag. The referenee magnitude for FDF in the B band is AIq = 5.48. The Seheehter LF for 
galaxies has been widely used to parametrise the LF in FDF as a funetion of the redshift. As an example 
in band B, a = —1.07 in the range 0.45 < 2 ; < 1.11 with —19.52 < M* < —18.80, see Figure 5 in [12]. 
We have maintained a = —1.07 but we make M* variable and speeify it in the eaptions of the figures. 
The distribution of the speetroseopie redshifts in the FDF is presented in Figure 6 and a eomparison 
should be made with the distribution of photometrie redshifts, see Figure 2 in [13]. 

Figure 7 presents the number of observed galaxies in the FDF eatalog at a random apparent magnitude 
and Figure 8 reports the theoretieal number of galaxies as funetion of redshift and apparent magnitude. 
Here we adopted the law of the rare events or Poisson distribution in whieh the varianee is equal to the 
mean, i.e. the error bar is given by the square root of the frequeney. An enlarged diseussion on the 
validity of this approximation ean be found in [35]. 

A eareful examination of Figure 3 in [13] gives the maximum frequeney of galaxies with well 
measured speetroseopie redshift in the FDF at z ~ 0.3. The mean redshift of galaxies as a funetion 
of the apparent magnitude for the FDF eatalog is presented in Figure 9, whieh shows an aeeeptable 
agreement between the data (empty stars) and the theoretieal values (full line). 

4.2. The zCOSMOS catalog 

The zCOSMOS bright redshift 10k eatalog, whieh eovers a solid angle Vt = 1.7deg^ or O = 
3.04617 lO^'^sr, eonsists of 9697 galaxies in the the interval [0,4] in redshift and range in the Ib band 
[15,23] mag, see [14]. The referenee magnitude for zCOSMOS in the Ib band is A4q = 4.08, see Table 
2.1 in [36]. The number of galaxies as a funetion of the redshift does not have a eontinuous behavior: 
rather, we are in the presenee of an alternating behavior of voids and relative maxima, see Figure 10. 
This nonhomogeneous spatial distribution of galaxies ean be made eontinuous by introdueing bigger 
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Figure 7. The galaxies of the PDF eatalog with 22.08 < m < 26.81 or 2.33jJ^ < 

f < 181.39;^^^ are organized in frequeneies versus speetroseopie redshift. The redshift 
covers the range [0,1.5] and the histogram’s interval is 0.18. The maximum frequency of 
observed galaxies is at z = 0.33, = 77.8, and the number of bins is 8. The full line is 

the theoretical curve generated by (z) as given by the application of the Schechter LF 
which is Eq. (48) with = 0.01 /Mpc^, M* = —17.78 and a = —1.07. 


photometric maximum 





Figure 8. The theoretical number of galaxies of the FDF catalog as afunction of redshift and 
apparent magnitude represented as a 3D surface, parameters as in Figure 7. 



Figure 9. Average observed redshift, (zobs), as a function of the apparent magnitude for 
the FDF catalog (empty stars) and theoretical full line for (z) as given by the numerical 
integration of Eq. (52). Theoretical parameters as in Eigure 7. 
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Figure 10. The galaxies of the zCOSMOS eatalog are organized in frequeneies versus 
speetroseopie redshift. The redshift eovers the range [0,1.2] and the histogram’s interval 
is 0.02. 


U) 



Figure 11. The galaxies of the zCOSMOS eatalog with 17.88 < m < 19.06 or 
803.43;^^^^ < f < 2392.36;^^^^ are organized in frequeneies versus speetroseopie redshift. 
The redshift eovers the range [0,1] and the interval in the histogram is 0.1. The error bar is 
given by the square root of the frequeney (Poisson distribution) . The maximum frequeney 
of observed galaxies is at z = 0.213, = 147.3, and the number of bins is 10. The full line 

is the theoretieal eurve generated by (z) as given by the applieation of the Seheehter 
LF, whieh is Eq. (48) with = 0.01 /Mpc^, M* = —20.88 and a = —1.07. 


intervals in the eomputation of the frequeneies, e.g., a histogram interval equal to 0.1. Figure 11 presents 
the number of observed galaxies in the zCOSMOS eatalog for a given apparent magnitude. 

The total number of galaxies, ean be eomputed with the integral represented by Eq. (51), see 
Figure 12. 

The mean redshift of galaxies as a funetion of the apparent magnitude for zCOSMOS is presented in 
Figure 13. 


5. The relativistic case 
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Figure 12. All the galaxies of the zCOSMOS eatalog, organized in frequeneies versus 
speetroscopie redshift. The redshift eovers the range [0,1] and the interval in the histogram 
is 0.1. The error bar is given by the square root of the frequeney (Poisson distribution) . The 
maximum frequeney of all observed galaxies is at 2 ; = 0.35, = 1864.65, and the number 

of bins is 10. The full line is the theoretieal eurve generated by as given by the 

numerieal integration of Eq. (51) with = 0.01 /Mpc^, M* = —18 and a = —1.07. 



Figure 13. Average observed redshift, (zobs), as a funetion of the apparent magnitude for 
the zCOSMOS eatalog (empty stars) and theoretical full line, (z), as given by the numerical 
integration of Eq. (52). Theoretical parameters as in Eigure 11. 
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The possibility of deriving an analytieal result for the number of galaxies as a funetion of the redshift 
in the relativistie ease is eonneeted with the availability of an analytieal expression for the luminosity 
distanee. We now use the same symbols as in [37] and we define the Hubble distance Dh as 

c 




We then introduee a first parameter 


M 


— 


Svr G po 

~YHf 


(53) 


(54) 


where G is the Newtonian gravitational eonstant and po is the mass density at the present time. A seeond 
parameter is 


Oa = 


Ac" 

Jiq ' 


(55) 


where A is the eosmologieal eonstant, see [31]. The two previous parameters are eonneeted with the 
eurvature by 

Om + Oa T — 1 


The eomoving distanee, Dq, is 


Dc = D 


H 


dz' 

W) 


where E{z) is the ‘Hubble function’ 

E(^z) = y/Oivi (1 + (1 + 

The transverse comoving distance Dm is 


D 


H 


sinh [y/n^ Dc/Dn] for 0^- > 0 


Dm= < Dc 

D. 


sm 


VW^\Dc/D 


An analytic expression for Dm can be obtained when Ha = 0: 


H 


for Dk = 0 
for Qk < 0 


Dm = Du 


2 [2 — Hm (1 — z) — (2 — Hm) a/I + Hm z 


(i- + 


for Ha = 0. 


A new form for Dm when Ha = 0 is 

Du sinh (2 arctanh 

Dm = -^- 


\/T—^ 


2 arctanh ( 

V y 


The luminosity distance is 


which in the case of Ha =0 becomes 


y/1 — Hm 

Dc = (1 + z) Dm 


Dc — (1 + 2 :) 


Du sinh ( 2 arctanh 


1 


\/l—O m 


2 arctanh 




M 


(56) 

(57) 

(58) 


(59) 


(60) 


(61) 


(62) 


(63) 
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and the distance modulus 


m — M = 25 + 5 

We now return to Dl = 
differential of the redshift is 


1 In ( 2 ^ ~ (1 ~ “ (2 ~ ‘*^m) V + l) 


ln(10)“^\“ 


(64) 


r: the ration between the differential of the luminosity distance and the 


dr = 


cf^M ^2 y/ f^M ^ “f 1 f?M — 2^ dz 


Hoy/Z + li?] 


(65) 


M 


This means that we have an analytical expression for the differential dr = f{z)dz when = 0. This 
analytical differential will be inserted later on in eqn.(67). The inverse relation between distance and 
redshift, now denoted by is 


zm = 


r T cVt 


M 


y/(2 rH q + c)f2M — 2 cflM + 2 y/(2 rH q + c) 


2 C^M 

The joint distribution in z and/for the number of galaxies in the relativistic case is 


dN 

dVtdzdf 


1 L L 


and its explicit value is 


where 


dN _ DNN 
dVLdzdf ~ DND 


DNN — 64 $ ^y/f2]vi + ~ 2 y/flM + 1 — 5" ^ 

T — 2 y/f2]vi + 1 — T 2^ \ 

I ( 2 ]y[ Zcrit I 

_ (-yflivi +in]y[+n]y[ -2 -ynjvi +1-0^+^)^ / . - y 

e (2 y/i^M+l + (2 m — 2j c^tt 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 


and 

DND = QuW^uz + lHo^L* . (72) 

Figure 14 presents the number of observed galaxies in the zCOSMOS catalog for a given apparent 
magnitude in the relativistic case; we adopted the value of (2 m = 0.286 because it is the concordance 
value, see [7]. 

The mean numerical redshift of galaxies as a function of the apparent magnitude for zCOSMOS is 
presented in Figure 15 for the relativistic case. 

The observed galaxy number over-density on cosmological scales up to second order in perturbation 
theory with all relativistic effects that arise by observing on the past lightcone are discussed in [38,39]. 
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Figure 14. The galaxies of the zCOSMOS catalog with the same parameters of Figure 11 are 
organized in frequencies versus spectroscopic redshift. The full line is the theoretical curve 
generated by {z) as given by the application of the Schechter LF in the relativistic case, 
which is Eq. (68) with = 0.01 /Mpc^, M* = —20.7, a = —1.07 and Dm = 0.286; 
=95.68 when the number of bin is 10. 



Figure 15. Average observed redshift, (zobs), as function of the apparent magnitude for 
the zCOSMOS catalog (empty stars) and theoretical full line, {z), as given by a numerical 
integration. Theoretical parameters as in Figure 14. 
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6. Evolutionary effects 

The main problem in modelling the LF as a funetion of the redshift is that the low luminosity galaxies 
progessively disappear. This observational faet ean be solved by adopting a truncated probability density 
function (PDF). The beta distribution is defined in [0,1] and the beta with scale PDF is defined in [0,6]. 
On introducing a truncation in the beta PDF at low values, we can model the observed LF as a function 
of the redshift, see formula (34) in [40]. Once the random variable X is substituted with the luminosity 
L, we obtain a new LF for galaxies, 4^, 

^{L)dL = {Lb-dL , (73) 

where 4^* is a normalization factor which defines the overall density of galaxies, a number per cubic 
Mpc. The constant Kl is 

^ (74) 

A = —a r (a + /9) 

B = Lb^-\V^{a, -/3 + 1; 1 + a; (a + (3) 

-L/+“-'r(i + a)r(/?) , 

and Lq, Lfo are the lower, upper values in luminosity and 2 Fi{a^h] c; z) is the regularized hypergeometric 
function [41,42]. The averaged luminosity, (L), is 

, (75) 

-dd 

where 


Aj<i — — Lb^ —[3 + 1,1 + 2 + a; ^)F (1 + + /?) 

+Lb“+^r (2 + a) r { 13 ) 

Bjy = (1 + a) r (1 + tt + /?) 

The relations connecting the absolute magnitude M, Ma and Mb of a galaxy to the respective luminosities 
are 

^ _ ^qOA{Mq-M) Lg ^ ^^0A{MQ-Ma) Lb _ ^Q0.4(MQ-Mi,) {16) 

Lq ’ Lq ’ Lq 

where Mq is the absolute magnitude of the sun in the considered band. The beta truncated LF in 
magnitude is 

^{M)dM = -Km'^*0A ^ X 

^^q- 0.4 Mj+0.4 Mq _ ^Q-0.4am+0.4MQ ^ 

10“0-4“"*+o.4Mq dM 


(77) 
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where 



(78) 


Am = aT{a + (3) 



-0.4M6+0.4Mq\^ ^ 


2Fi(cr, —(3 + 1; 1 + cr; 


^Q—0.4Ma+0.4MQ 


^g— O.4M(,+O.4M0 



-0.4Ma+0.4M. 


'0)“r(a + « 



—O.4Af^H-O.4Af0 


/3+q;—1 


r(i + a)r(/?) 


This LF eontains the five parameters a, (3, Ma, Mb, \I/* whieh ean be derived from the operation of 
fitting the observational data and Mq whieh eharaeterize the eonsidered band, see [43] . The number of 
variables ean be reduced to three once Ma and Mb are identified with the maximum and the minimum 


absolute magnitude of the considered sample. A further reduction of parameters can be realized in the 


case of a well defined catalog of galaxies, e.g., zCOSMOS, where Ma=-23A1 mag, a = 0.01, (3 = 5 
and Mq=4.08. The low luminosity bound (high magnitude) can be modelled in the classic case by 
extracting the absolute magnitude from Eq. (28) which represents the distance modulus for tired light by 



(79) 


where mi is the limiting apparent magnitude, which for zCOSMOS is rriL = 23.2 mag. With the 
above choice of parameters, the observed LF for zCOSMOS as a function of the redshift has only one 
free parameter, 4/*, which can be easily derived from the fit of the histograms. The observed LF for 
zCOSMOS can be built by adopting the following algorithm. 

1. A value for the redshift is fixed, z, as well as the thickness of the layer, A 2 ;. 

2. All the galaxies comprised between z and Az are selected. 

3. The absolute magnitude can be computed from Eq. (28) which represents the distance modulus for 


tired light. 


4. The distribution in magnitude is organized in frequencies versus absolute magnitude. 


5. The frequencies are divided by the volume, which is L = Ovrr^Ar, where r is the considered 
radius, Ar is the thickness of the radius, and O is the solid angle of ZCOSMOS. 


6. The error in the observed LL is obtained as the square root of the frequencies divided by the 
volume. 


Ligures 16, 17, and 18 present the LL of zCOOSMOS as well the fit with the truncated beta LL at 
z = 0.2, 2 : = 0.5, andz=0.1, respectively. 
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Figure 16. The luminosity function data of zCOSMOS are represented with error bars. The 
continuous line fit represents our beta LF (77), the parameters are z = 0.2, Az=0.05 and 
NDIV =8, which means = 5.35. 



Figure 17. The luminosity function data of zCOSMOS are represented with error bars. The 
continuous line fit represents our beta LF (77), the parameters are z = 0.5, Az=0.05, and 
NDIV =8, which means =16.71. 



Figure 18. The luminosity function data of zCOSMOS are represented with error bars. The 
continuous line fit represents our beta LF (77), the parameters are z = 0.7, Az=0.03, and 
NDIV =4, which means =8.76. 
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Figure 19. The luminosity function data of zCOSMOS are represented with error bars. The 
continuous line fit represents our beta LF (77) in the relativistic case. The input parameters 
are = 0, = 0.286, z = 0.7, Az=0.03 and NDIV =4, which means =8.78. 


In the relativistic case, we can extract the absolute magnitude from Eq. (64), which represents the 
distance modulus when (Ia = 0: 


Ml, = mi — 5 


In (2) 

In (10) ^ln(lO) 


-5 


In 


C (2 — [1 — z) — (2 — I^m) a/ zf2M + l) 


Hna 


-25 


(80) 




Figure 19 presents the LF of zCOOSMOS as well as the fit with the truncated beta LF when = 0.7 in 
the relativistic case. 


7. Conclusions 

Results: 

A nonlinear formulation of Hubble’s law allows determinaing an old/new relation for the distance as 
a function of the redshift, see Equation (7). This distance, when inserted in the definition of the joint 
distribution in z and/for the number of galaxies, allows the determination of a new N — z relation, see 
Eq. (48). Two photometric tests were done for 263 galaxies belonging to the FDF catalog and for 9697 
galaxies belonging to the zCOSMOS catalog. The first test is dedicated to the photometric maximum in 
redshift, for which is possible to derive an analytical expression as a function of the flux, see Eq. (49), or 
the apparent magnitude, see Eq. (49); this results in Figures 7 and 11. 

A second test is dedicated to the average redshift, for which a numerical integration of Eq. (52) should 
be done; it results in Figures 9 and 13. 

The same formalism can also be applied to the relativistic case once the relativistic luminosity distance 
is given, see Eq. (63). In the standard cosmology or relativistic case the the joint distribution in z and 
/for the number of galaxies is given by Eq. (68). A comparison between the Euclidean and relativistic 
model can be made on the photometric maximum as represented by Figures 11 and 14. The test gives 

= 147.3 for the Euclidean case, as represented by Eq. (47) and = 94.27 for the relativistic case as 
represented by Eq. (68), but large oscillations are present in the observed frequencies and therefore the 
definitive answer is remanded to future efforts. 

The observed LF for galaxies can be modelled by a truncated beta LF, see Eq. (77). This new LF with 
an appropriate choice of parameters has only one free parameter, which is the number of galaxies per 
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cubic Mpc, and this parameter deereases with the redshift. The high magnitude bound, M^, ean be 
modelled both by a Euclidean model as given by Eq. (79) and by a relativistie model, = 0, as given 
by Eq. (80). As an example, the third test for the observed EE for galaxies dX z = 0.7 gives =8.76 for 
the Euelidean ease and =8.78 for the relativistie ease when I^a = 0 and = 0.286. 

Generalizated tired light: 

The presenee of the factor for adjustable tired light, see Eqs. (29) and (30), poses the problem of 
its determination. Eq. (30), whieh represents the distanee modulus, ean be ealibrated on the database of 
supernova (SN) of type la. A eareful determination of and Hq ean provide a better determination of 
the Malmquist bias, as represented by Eigures 3 and 4, whieh present a laek of galaxies just above the 
red lines. 

Tired light versus GR 

The new distance modulus as represented by Eq. (30) requires a eareful eomparison with the standard 
ECDM. In the ease of ECDM with (Ia = 0, an analytieal solution for the luminosity distanee exists and 
allows the determination of the relation between the differential of the distanee and the differential of the 
redshift, see Eq. (65). 

In the ease of IIa 0, an analytieal expression for the luminosity distanee does not exist, and as a 
eonsequenee we do not have at the moment of writing a relation between the differential of the distanee 
and the differential of the redshift. The use of the Fade approximant ean produee analytical results for 
the luminosity distanee, see [44,45]; this approximation will be the subjeet of future researeh. 

The cells 

The eellular strueture of our universe, see as an example Eigure 10, is now the greatest ineonvenienee 
to the application of the eontinuous models for the number of galaxies as a function of the redshift, and 
perhaps an explanation for why the theoretical lines do not fit the data, as an example see Eigure 12. We 
briefly reeall that there is an actual debate on the dimension of the universe which is modelled by 77 oc 
R^, where N is the number of galaxies, R the radius of the eonsidered sphere, and D the dimension. An 
homogeneous universe means 79 = 3. In the eoneordanee model, D makes a transition to D=3 at seales 
between 40 and 100 Mpe, see [46]. An aceurate analysis of 2MASS Photometric Redshift catalogue 
(2MPZ), shows an agreement with the standard eosmologieal model; the homogeneous regime is reaehed 
faster than a elass of fractal models with D < 2.75, see [47]. As an example of non-homogeneity the value 
79 = 1.87 has been reported in [48] where the 2MASS Redshift Survey eatalog was analysed. 
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